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Time- and State-Dependent Input Delay- 
Compensated Bang-Bang Control of a 
Screw Extruder for 3D Printing 

Mamadou Diagne^, Nikolaos Bekiaris-Liberis^, and Miroslav Krstic^ 


Abstract —In this paper a delay-compensated Bang-Bang con¬ 
trol design methodology for the control of the nozzle output flow 
rate of screw-extruder-hased 3D printing processes is developed. 
The presented application has a great potential to move heyond 
the most commonly used processes such as Fused Deposition 
Modeling (FDM) and Syringe Based Extrusion (SBE), improving 
the huild speed and the 3D parts accuracy. A geometrical 
decomposition of the screw extruder in a partially and a fully 
filled regions (FEZ and FEZ) allows to describe the material 
convection in the extruder chamber by a ID hyperbolic Partial 
Differential Equation (PDE) conpled with an Ordinary Differen¬ 
tial Eqnation (ODE). After solving the hyperbolic PDE by the 
Method of Characteristics (MC), the coupled PDE-ODE’s system 
is transformed into a nonlinear state-dependent input delay 
system. The aforementioned delay system is extended to the non- 
isothermal case with the consideration of periodic fluctuations 
acting on the material’s convection speed, which represent the 
effect of viscosity variations due to temperature changes in the 
extruder chamber, resulting to a nonlinear system with an input 
delay that simultaneously depends on the state and the time 
variable. Global Exponential Stability (GES) of the nonlinear 
delay-free plant is established under a piecewise exponential 
feedback controller that is designed. By combining the nom¬ 
inal, piecewise exponential feedback controller with nonlinear 
predictor feedback the compensation of the time- and state- 
dependent input delay of the extruder model is achieved. Global 
Asymptotic StablUty (GAS) of the closed-loop system under 
the Bang-Bang predictor feedback control law is established 
when certain conditions, which are easy to verify, related to 
the extruder design and the material properties, as well as to 
the magnitude and frequency of the materials transport speed 
variations, are satisfied. Several simulations results are presented 
to illnstrate the effectiveness of the proposed control design. 


I. Introduction 

Additive Manufacturing (AM) has a promising future and 
demonstrates its effectiveness in various applications involving 
tissue engineering ||T|, Q, chemical engineering Q, thermo¬ 
plastics Q, metal Q and ceramic @ materiars fabrication. 
Functional 3D objects with complex geometrical shape can be 
produced in a short time without the needs of tools thanks to 
the Computer Aid Design (CAD) that drastically reduces the 
products development procedure. Currently, the most popular 
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plastics 3D printers are based on FDM Q-Q and SBE Q, 
p0| technologies (Fig[T]i. 

In these processes, biodegradable polymers are transported, 
heated and pressurized in an extruder chamber before being 
dropped on a platform, one horizontal thin layer at a time, until 
the complete 3D part is built such that it closely resembles 
the original CAD model. One of the crucial point that is not 
commonly addressed in the existing literature of extrusion- 
based 3D printing is controlling of the start and stop of 
extrusion-on-demand. A hybrid extrusion force-velocity mod¬ 
eling and tracking control for the fabrication of functionally 
graded material parts is developed in GD using a first order 
differential equation that describes the plunger dynamic in a 
SBE process. Some extents of that approach are proposed by 
1121 with a robust tracking of the extrusion force to recover 
constant flow disturbances whereas GD considers an unknown 
transfer function gain with an adaptive control strategy. Several 
issues regarding on EDM are discussed in GD, GH and 
references therein, including the potential clogging due to 
agglomerate formation at the nozzle, appearance of bubbles, 
density inhomogeneity, tracking of short time-scale process 
variations, and prediction of anomalies such as material over¬ 
flow and underflow for diverse applications. Thermal control 
is left out of most prior studies which are essentially based on 
empirical models. 



Fig. 1. (a) FDM, (b) SBE, (c) SE processes 0. 

In this paper, we are interested in the flow control issues 
related to the recent advances of 3D printing technology for 
which a Screw Extrusion (SE) process is utilized. In SE, the 
rotating screw allows a continuous feeding mechanism and 
generates a sufficiently high pressure in the extruder chamber, 
increasing, as a result, the printing speed. In addition, the 
screw motion extends the mixing capabilities of the system, 
and thereby, reduces drastically the risk of potential clogging 
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at the nozzle while improving the homogeneity of the extruded 
filament 0^ HD- The SE process with granular material 
moves beyond the restrictions of EDM and does not require 
hlament-shaped raw materials to operate. Consequently, it 
enables the processing of a broader range of raw materials and 
permits an easy recycling of wasted plastic during extrusion 
0, m- In such processes, the need to control the start and 
stop of the extrusion process on demand calls for advanced 
control methodologies that are capable of enhancing the final 
product’s quality in an industrial level. Even if experimental 
results demonstrate the effectiveness of SE d), 0, the 
challenging control problems arising in such applications are 
actually poorly investigated. 

In the present article, a generic and dynamical model of 
a homogeneous melt SE process derived from mass and 
momentum balance laws p8) is used for the design of 
a delay-compensated Bang-Bang controller which permits a 
fast and accurate control of the flow at the nozzle output. The 
model consists of a ID Partial Differential Equation (PDE) 
that is defined on a time-varying spatial domain whose dy¬ 
namics obey to an Ordinary Differential Equation (ODE). The 
transformation of the coupled PDE-ODE system into a state- 
dependent input delay system, which describes the dynamics 
of the material convection in the extruder chamber, is achieved 
after solving the PDE by the Method of Characteristics (MC) 
d), |Tg. In order to also account for potential periodic 
fluctuations of the materials transport speed when processing 
granular pellets d), due to the thermal energy that is supplied 
into the system from the heater of the extruder and due to the 
mechanical shearing effect by the rotation of the screw, the 
state-dependent input delay model is extended to a nonlinear 
system with an input delay that depends simultaneously on 
the state and the time variable (see 1^ , | |2^ and | |22) for the 
treatment of systems with time- and state-dependent delays). 

In p9) , a delay-compensated Bang-Bang control law is 
developed for the control of the nozzle output flow rate of 
an isothermal screw extrusion process, achieving GES of the 
delay-free plant at any given setpoint. By combining the 
nominal, piecewise exponential feedback controller p9) with 
nonlinear predictor feedback, which is extended from the 
state-dependent input delay case |20| to the case in which 
the vector held and the delay function depend explicitly on 
time, the compensation of the time- and state-dependent input 
delay of the non-isothermal screw extrusion model is achieved. 
GAS of the closed-loop system under the delay-compensated 
Bang-Bang controller is established when certain conditions, 
related to the extruder design and the material properties, 
as well as to the periodic Huctuations, are satished. Several 
simulations results are presented including the case in which 
there is uncertainty in the value of the periodic variations of 
the material’s transport speed. 

This paper is organized as follows: The screw extruder 
mechanisms and the bi-zone model of the extruder consisting 
of the transport PDE coupled with the ODE for the moving 
interface is discussed in Section [II] In Section [111] the transfor¬ 
mation of the coupled PDE-ODE system into a state-dependent 
input delay system by computing the PDEs solution by the 
MC is presented and it is then extended to a nonlinear system 


with a time- and state-dependent input delay. The control of 
the delay-free plant with a piecewise exponential Bang-Bang- 
like control law is described in Section |IV] In Section |V] 
we design the predictor feedback control law for nonlinear 
systems with time- and state-dependent delay acting on the 
input. The application of the predictor feedback control law to 
the screw extruder model is presented in Section VI The paper 
ends with simulations, including a discussion on the robustness 
properties of a state-dependent input delay compensator to 
time- varying perturbations acting on the vector field and the 
delay function, in Section [Vn| 


II. 3D Printing Based on Single-Screw Extruders 

A. Extrusion process description and structural decomposition 
of the extruder into a partially and a fully filled zone 

A screw extruder is divided into one or several conveying 
zones (transport zones), melting zones (for material fusion) 
and mixing zones in which the extruded melt is submitted to 
high pressure, before its eviction through the nozzle |T3,|Tg, 
p4|-p0|. The net flow rate at the extruder nozzle depends 
mainly on the material flow in the longitidunal direction 
given by ID heat and mass transport equations (3T), (32). 
Another particularity of these processes is that they can be 
divided in geometric regions which are partially and fully filled 
called PEZ and EEZ, respectively (Eig. |^. The PEZ which is 
submitted to an atmospheric pressure is a conveying region 
and the flow in the EEZ is determined by the pressure gradient 
building-up in that region due to the nozzle resistance. These 
two zones are coupled by an interface which moves according 
to the volume of material accumulated in the EEZ. Basically, 
the moving interface is located at the point where the pressure 
gradient passes from zero to a non null value. 



A 


L 


Eig. 2. Bi-zone model of a screw extruder. 


B. Mass and momentum balance of an extrusion process 
1) Mass balance of the PFZ : The PEZ is dehned on the 
time-varying spatial interval {x{f),L), x{f), being the length 
of the EEZ and L the extruder length (Eig. ig. The mass 
balance in this area can be expressed using the fraction of the 
effective volume between a screw element and the barrel (14ff) 
which is occupied by the extruded material, namely, the hlling 
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ratio u. Considering an incompressible homogeneous mixture 
with constant density po and viscosity p, the following mass 
conservation equation is deduced 

dtu{z,t) = e (K+, (a;(f),L)) (1) 

u{L,t) = U{t), (2) 


where Nq is the constant screw speed and ^ the uniform 
pitch of the screw. The boundary condition u{L, t) is defined 
assuming the continuity of the flow at the inlet {z = L} 


U{t) = 


FUt) 

PoNoVeB ’ 


(3) 


where Fi„(f) is the feeding rate. Physically, the term pqNqVsB 
in 0 is the maximum pumping capacity of the screw. 

2 ) Momentum balance of the FFZ: The FFZ whose filling 
ratio is equal to one is defined on the spatial domain (0, x{t)), 
where the coordinate {z = 0} is the extruder’s end. The FFZ 
flow depends on the pressure gradient that appears in this 
region, resulting to backward or forward flow. The momentum 
balance which is derived from Navier-Stokes equations under 
stationary conditions yields the pressure gradient 


Bpo 


(4) 


for all {t,z) € (R~^, (0,x(t))), where B is a coefficient 
of pressure flow. The net flow rate Fd{f), in the case of 
a Poiseuille flow is expressed with the help of the nozzle 
conductance Kd, the viscosity rj, and the pressure at the nozzle 
P{0,t) as 


Fdit) = f AP(<), 

AP{t) =P{0,t)-Po. 


3) Mass balance of the FFZ : The FFZ mass balance leads 
to an ODE which describes the time evolution of its length. 
This length denoted by x{t) determines the location of the 
small transfer region that is assimilated to the point at which 
the pressure changes from the atmospheric pressure Pg to a 
different value |T8|, 113, mg, 130), fn) 

dxit) _ po^NoVesu(x{t),t) - Fdit) 
dt poSesi^ - u{x{t),t)) 

where 5'eff is the available section and Ves = 

4) Coupling relations at the PFZ-FFZ interface : The cou¬ 
pling condition is formulated imposing the pressure continuity 
at the spatial coordinate x{t) 


P{x-,f) = P{x+,t) = P^. (7) 


Integrating the pressure gradient equation (|4|, the net flow rate 
defined in (|5]l is written as 


p ^ KdVeffNoPox{t) 
^ Bpo + Kdx{t) 


( 8 ) 


Substituting (|5]l in (|g and using ([8]), equation 0 for the length 
of the FFZ is written as 


dx{t) 


Kdx{f) - {Bpo + Kdx{t)) u{x{t), f) 
{BpqF Kdx{t)){l-u{x{f),t)) 


III. From Mass Balance Equations of the 
Extruder to a Delay System 


A. Isothermal delay system model 

The bi-zone model ([T]i, (|2]l, and (|9]l can be reduced to 
a nonlinear state dependent-input delay system p0[ . The 
characteristic solutions of ([T]i with respect to the boundary 
condition (|2| are 

u{z,t)=u[t-^y ( 10 ) 

Substituting ( [TOl i into 0, we derive the following nonlinear 
system 


x(f) =^No 


Kdx(t) 


{Bpo + Kdx{f)) {l-U {t- Ds{x(t)))) 
U{t- Dsixlf))) 




P(f)G[0,l) (11) 


The state-dependent input delay function is denoted as 

Dsixit)) = (12) 

A detailed derivation of the ODE 0 and an extensive 
description of the screw-extruder model for 3D printing is 
given in ng. 


B. Delay system representation in a non-isothermal case 

In this section, we propose the following extension of the 
state-dependent delay model o to account for the viscosity 
changes due to the temperature variations occurring in the 
physical system when processing plastic pellets 



[ {l + e2x{t))(l-U( 
U{t-D{x{f))) 

t-D{t,x{f)))) 

, (/(<)£ [0,1), (13) 

{l-U{t-D{fx{t))))_ 


where 


and. 


01 

02 


^No, 

Kd 

Bpo' 


(14) 

(15) 


D{t,x{t)) = 


L — x{f) 


c{f) 

c{t) = 0i(l -f ecos{ujt)), 


(16) 

(17) 


where, e < 1 is a positive constant and w is the mean value of 
the angular frequency of the periodic fluctuations. Our choice 
for the non-isothermal model is motivated by the fact that 
the expansion of granular material into a plastic state due to the 
thermal effect leads implicitly to periodic fluctuations of the 
convection speed, namely, 0i due to the viscosity variations 
m- Moreover, some nozzle “instabilities” phenomena may 
appear as short periodic distortions of the extrudate, due to 
the viscoelastic properties of the fluid, with magnitude smaller 
than one HD- 


dt 
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IV. Control of the Delay-Free System with a 
“Bang-Bang” Control Law 


A. Open-loop stability 


The starting point of the delay system controller design 
consists of the construction of a nonlinear control law that 
stabilizes the delay-free system 


i{t) =c{t) 


02x{t) 


{1 02x{t)) {1 - U (t)) 

Uit) 




t/(f)e[o,i). (18) 


For e < 1, the time-varying speed of the material transport 
c(f) is strictly positive and the open-loop stabilizing control 
law of the delay-free plant (18 1 is given by 

e2x* 


j{x*) = 


1 -f 92X* 


'ix* G [0,L), 


(19) 


for the physical parameters of the extruder satisfying 02 L < 1. 
This statement is directly derived considering the Lyapunov 
function V = |e(f)|, where e{t) = x{t) — x*. 


B. “Bang-Bang" controller design with piecewise exponential 
functions 

For the feedback stabilization of ( |T8] l, we consider two 
exponential functions 

• for x{t) < X*, a left-exponential function 

vfx,x*) = ’ ^^0) 


where ai{x*) > 0 is the gain of the left exponential con¬ 
trol law. The function ( |20| ) takes values in Umax], 

where Umax < 1 is the maximal value of the inlet 
filling ratio, namely the maximal feeding capacity of the 
extruder. Therefore, ui(0) = Umax, allows to set the inlet 
flow at its maximum capacity for a rapid refill action 
when the extruder is empty. 

• for x{t) > X*, a right-exponential function 


Ur(a::, x*) = v{x*) 


v{x*) 


2 _ g-ar(a:*)(a;-x*) 

2 _ Q—a-r{x*){L—x*) ’ 


( 21 ) 


where ar{x*) > 0 is the gain of the right exponential 
control law. The function (|2T| belongs into the interval 
[0, v{x*)] and the control action stops radically the flow 
when the extruder is completely filled, namely, Ur(L) = 

0 . 


controllers and respectively into the differentiable 
piecewise exponential feedback law ( [2^ . The slope function is 
defined as S{x) = ^ (the minus sign is conventional). 

More precisely, the key point of the design is to define a 
free parameter that may be specified by the user as the 
value of slope function at the equilibrium S{x*), under some 
restrictions that will be emphasize in this section. It is clear 
that, equiting the assigned value S{x*) to both left and right 
slope functions of ( |20) l and we can easily derive the 

following relations 


Six*) 

Six*) 


aiix*)iVm.ax - Vjx*)) 

^ _ ^—ai{x*)x* 

arix*)vix*) 

^ ^ — ar{x*){L — x*) 


(23) 

(24) 


The equations ( |2^ and are both transcendental and 

admit numerical solutions namely the suitable exponential 
parameters needed to the left and to the right of the setpoint 
to achieve the differentiability of the controller (22 1 . These 
solutions ai(x*) > 0 (respectively, a^ix*) > 0) exist if 
the linear and exponential functions of ai(a:*) (respectively, 
ar(a;*)) have a strictly positive intersection. Consequently, the 
desired slope function Six*) should be above some minimum 
value denoted S^ninix*) , for any given equilibrium in the 
physical domain (0, L). More precisely, the two equations in 
( |2^ have strictly positive solutions if at the origin (oi = 0 and 
Or = 0), the slope of their linear part is less than the slope of 
their exponential part respectively. 


For the left exponential slope function (23 i, we define the 
linear and the exponential functions 


Itpiiaiix*)) = aiix*)ivma^ - vix*)) 

a solution of ( [2^ should satisfy 

j'0i(ai) =Mai), 

t dv>i( 0 ) , d<i>i( 0 ) 

da] da] 


(25) 


(26) 


It follows that for all x G (0,a;*), Six*) should satisfiy 
the inequality 


Six*) > 


Vmax - vix*) 


(27) 


Decomposing the right exponential slope function (24i 
into 


C. Extension of the “Bang-Bang” control law on the whole 
domain (0, L) 

Next, we introduce the characteristic function of the do¬ 
mains [0, X*] and [x*, L] and write the extended control law 
as 

vix, X*) = viix, x*)hix* — x) -\- v-cix, x*)hix — x*), (22) 

where h is the Heaviside function. 

A continous slope function at the setpoint x* denoted by 
Six*) is imposed to extend the left and the right exponential 


I Ipriarix*)) = arix*)vix*) 

we deduce that a solution of ( [24| should satisfy; 

itprittrix*)) = friar ix*)), 

t dtllrjO) , d<t>r jo) 

\ dflr dar 

Hence, for all x G (x*,L) 


5'(x*) > 


vix*) 

L-x* 


(28) 


(29) 


(30) 
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Finally, the minimal value of the setpoint slope 5'inin(a;*) 
above which the gains ai(a;*) and a^ix*) ensure the differen¬ 
tiability of the extended control law ( 22 i on ( 0 , L) is given 
by 


Su 


[X = 


^ + 1 


max • 


1 


1 


- 1 



The speed of the controller or its “agressivness” increases with 
the rise of the setpoint slope S{x*). As it is illustrated in Fig. 

with the characteristics of the control law for the setpoints 
X* = 0.02m and x* = 0.16m with different values of the 
setpoint slope value S{x*). 



Hence, 

V = e{t) sgn (e(f)), 


and with the help of (36i, by choosing U(t) 
is written as 

V = -a{t)V - 


(38) 
v{x(t), X*), 

(39) 


where a{t) > 0 and /3(f) > 0 for all x G [0, L) and U G [0,1). 
The functions a{t) and /3(f) are given by 


a{t) 

I3(t) 


c{t)02 

1 + d 2 {e{t) + x*) 

1 

(1 - v{e{t) + X* ,x*)) 
c(f) |'u(e(f) + X*^x*) — 'u(a;*)| 
(1 - v{x*)) (1 + 6 > 2 (e(f) + X*)) ’ 


(40) 


(41) 


where, c(f) > (1 — e), for all f > 0. Therefore 




1 - 6»2L 


(42) 


From ( [42| the closed-loop system is exponentially stable at 
X* G [0, L) for all ccq G [0,L). ■ 


V. Predictor Feedback Control eor Nonlinear 
Systems with Time- and State-Dependent Input 
Delay 


Fig. 3. Control function for different Smin setpoints. 


Theorem 1, For any setpoint a;* G [0, L) and for any choosen 
setpoint slope S{x*) G K satisfying S{x*) > S'i„in(a;*), 
where S'min(a^*) A given by {31), taking the control gains 
(ai{x*), ar{x*)) as solutions of 

ai(:T*)(t’max - v{x*)) - 5(x*)(l - = Q, (32) 

a,(a;*)v(a:*) - 5(x*)(l - = 0. (33) 

The closed-loop system consisting of with an initial 
condition xq G [0, L) and the extended control law @-(|22]) 
is GES at X = x*. 

Proof: We rewrite the delay-free plant ( [T8| ) as 

c(f) (C/(f)-t;(x*)) 


e(f) = 


(1 - v{x*)) (1 -I- 6 » 2 (e(f) -I- a;*)) (1 - [/(f)) 
c(f) 02 e(f) 


(34) 


1 + 6 ' 2 (e(f) -I- X*) ’ 
where, e(f) = x{t) — x*. The control law ( [2^ is written as 

v{x(t),x*) = v{e{t)X*,x*). (35) 


The extended control law ( [22| l is a decreasing function of x 
and consequently for all x{t) G [0, L), 

sgn (u(e(f) -f x*,x*) — v{x*,x*)) = — sgn (e(f)). (36) 

Moreover, 0 < v{x,x*) < Umax < 1 and v{x*,x*) = v{x*) 

is defined as the setpoint open-loop control Next, we 
introduce the following Lyapunov function 


V=\eit)\. 


(37) 


A. Predictor feedback design 

We consider the following nonlinear system with a time- 
and state-dependent input delay 

= f (43) 

0(f) = t - D{t,x{t)), (44) 

where x G K", U : [tg — D (fg, x(fo)), oo) —M, f > fg > 0, 
D G (K+ X IR";K_|_), and / : i?_|_ x K” x K —>• K” is 
locally Lipschitz with f {t, 0,0) = 0 for all f > 0 and there 
exists a class /Coo function a such that 


|/(f,a:,[/) I < d(|a;|-f |[/|). (45) 

The predictor feedback control law for system ( |4^ , ( |44l i is 
[/(f) = K{(j{t),P{t)), (46) 

where, for all f — D{t, x{t)) < 0 <t 

/(a(s),P(s),[/(s)) 


P{0) = x{t) + J^ 

a{0) = ^ ^ 


t-D(t,x(t)) 1 - P{(^{s),P{s),U{s)f^' 
1 


t-D{t,x(t)) F{a{s),P(s),U{s)) 


(47) 
ds, (48) 


and 


Fia{0),Pi0),U{0)) ia{0),P{0)) 

+ ^{ai0),Pi0)) 

X f{a{0),P{0),U{0)). (49) 

When simulating the predictor feedback controller ([47]i 


(46 1 , at each time step the ODE for the system (43 1-(|44|) 
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must be solved (using, for example, a simple Euler scheme) 
and the length of the delay must be computed (for example 
as the integer part of N(i) = say -/V(i), where t 

is the discretization step). The predictor is then computed by 
integrating simultateously the two integral relations ( [47| ) and 
( |48| l at each time step, using a numerical integration scheme. 
For instance, with the left endpoint rule of integration we get 


k—i—1 

P(z)=a;(z) + r ^ 

k^i-N{i) 
k—i — 1 

a{i)=i + T 

k=i-N{i) 


f{a{k),P{k),U{k)) 

l-F{a{k),P{k),U{k)y 

1 

{l-F{a{k),P{k),U{k)y 


(50) 

(51) 


The prediction of the state at the time when the current 
control will have an effect on the state is defined as 


Assumption guarantees that system x = is 

strongly forward complete with respect to w. 

Assumption 2. There exist a locally Lipschitz function k € 
([fojOo) X ]R";IR) and a function p £ /Coo such that the plant 
X = f (t, x{t), K(t, x{t)) +uj(t)) is input to state stable with 
respect to uj and k is uniformly bounded with respect to its 
first argument, that is, 

\K{t,x)\ < p{\x\) V t>to. (59) 

Assumption 3. D G C^{R+ x ^ and ^ are 

locally Lipschitz (to guarantee the uniqueness of solutions), 
and there exist class /Coo functions p 4 , p^, pg and non¬ 
negative constants ci, C 2 , C 3 , with C 3 < c, for some 0 < c < 1 , 
such that 


P(t)=x{t + D{<j{t),P{t))), (52) 

where the prediction time is defined as 

a{t) = t-\-D{a{t),P{t)), (53) 


which is derived from the inversion of the time variable 
t -)■ t- D{t,X{t)) + D{a(t),P{t)), f^, (2^. 

Differentiating ( [52] i, ( [5^ and using ( |4^ we arrive at 

= /(„((),, 54 , 

and 

= l-F(a{t),Pit),U(t)) 

where F is defined in ( |49] ). Finally, the implicit integral 
relations •ED and ([48) are derived by integrating and 
on the delay interval 

The key point of the predictor feedback design is the 
feasibility condition defined as 




-^{a{e),P{9)) 

+ ^ (0(9), P{9)) f{ai9), P{9), U{9)) < c,(56) 


for all 9 > to — D{tQ,x(tQ)) and some c € (0,1). Condition 
( |56| l guarantees that the feedback control action can reach the 
plant, namely, the delay rate is bounded by unity, and that 
the denominator of the predictor ( |47| i and the prediction time 
( |48] l is positive. We refer the reader to p0| for details on 
the predictor feedback control design and analysis for systems 
with state-dependent input delay. 


B. Stability analysis 

Assumption 1. There exist a smooth positive definite function 
R and class /Coo functions pi, p 2 and p^ such that for the 
plant X = f{t,x,w), the following hold 


D{t,x{t)) 

dD 

dD 


< 

< 

< 


Cl + P4{\x\), 

C2 + P6i\x\), 

C3+ d5i \ y )- 


(60) 

(61) 

(62) 


The definitions of strong forward completeness and input- 
to-state stability are those from p^ , and p5] , respectively. 

Theorem 2. Consider the closed-loop system consisting of the 
plant (|S i and the control law ( |47P -( |46| ). Under Assumptions 
m anc/|^ there exist a class K. function V'RoA ^ class 
ICC function j3s such that for all initial conditions for which 
U is locally Lipschitz on the interval [fg — D{to, x{to)),to) 
and which satisfy 


D{to) < V'RoaCc- C 3 ), (63) 


for some 0 < c < 1 , where 

n(/) = |a;(/)| -I- sup \U(9)\, (64) 

t — D{t,x{t))<9<t 

there exists a unique solution to the closed-loop system with x 
Lipschitz on [tg, cxd), U Lipschitz on (/g, cx)), and the following 
holds 


Cl{t) < l3s{Cl{to),t-to), (65) 

for all t > /g. Furthermore, there exists a positive constant 7 
such that for all t > /g. 


D{t,x{t)) < 7 

dD{t,x{t)) 

-;- < c 

dt 


( 66 ) 

(67) 


di{\y) < < T2 {\x\) (57) 

dR(t,x) dR(t,x) X X /I IX 

— dt --b^iadwl), (58) 

for all (x, Lo)'^ G and t > fg. 


Proof of Theorem Estimates ( |65] l, ( | 6 ^ , and ( |67] i follow 
by directly applying Lemmas 1-8 from | |20) (see the Ap¬ 
pendix). Existence and uniqueness of a solution x Lipschitz 
on [0,oo) follows from the proof of Theorem 1 in (page 
7). It remains to show that U is Lipschitz on (/g,(X)). Since 
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U{t) = K (cr(i), P(t)) and 

f {(T{t), P{t), K{(T{t), P(t))) 


P{t) = 
a{t) = 


1 - F{a{t),P{t),K {a{t),P{t))) ’ 
1 

1 - F{a{t),P\t),K (a{t),P{t)))" 
F{a{t),P{t))=^D{a{t),Pit)) 

X / {a{t), P{t), K P{t)))) 
iait),Pit)), 


+ 


at 


( 68 ) 

(69) 


(70) 


for t > to, the Lipschitzness of k and /, and (|56|l 


iitz 


ensure that the right hand-side of ( |68| ) and ( |69] l are Lipsc 
and consequently (P, cr) £ {C^{to,oo) x C^(fo)Oo)). From 
the Lipschitzness of k, it follows that U is Lipschitz. ■ 


VI. Application to the Extrusion Process Model 

From now, we recall the predictor feedback ([47|-([46| for 
the compensation of the time- and state-dependent input delay 
in system ([T^ that we rewrite formally as 

i{t) = f {t-D{t,x{t)))) ,(7l) 

f {t,x{t),U{t)) = -c{t)r{x{t),U{t)), (72) 


where D{t,x{t)) and c{t) are dehned in (16i and (HI, 
respectively, and 

02Xit) U (t)) 


T{x{t),Uit)) = 
Uit)e [0,1). 


il + 92xit))il-U{t)) (l-Uit)) 


(73) 


The predictive feedback controller based on the piecewise 


exponential feedback law (221 is given by 


U{t) = V (Pit)), (74) 

f{a{s),P{s),Uis)) 


P{e)=xit) + J^ 


t-D(t,x{t)) F{a{s),P[s),U{s)) 


•ds 


a{e)=t + J' 


(75) 
-ds, (76) 


t-D(t,x(t)) 1 - P(ct(s),P(s),17(s)) 

for all t — D{t,x{t)) < 0 < t. The function F dehned in 
(49 I for the system (|7T]i-(|72|i is computed with the help of 
(16 i-(fT7]) as 

Oieuj sin(wcr(f))(L — P{t)) 


F{a{t),P{t),U{t))=- 


c^it) 

T{P{t),U{t)), 


(77) 


where r{P{t), U{t)) is dehned in ( |73l ). 

The parameters ai{x*) and ar{x*) of the feedback control 
law are the solutions of ( [23] l and ( [24| l for an assigned 
slope function value at the set point that satishes (|^. P{t) = 
X {t + D{(j{t), P{t))) is the prediction of the state at the time 
when the current control will have an effect on the state. 
Recall that the implicit integral relation ( |75l l is derived from 
the inversion of the time variable t ^ t — D{t,x{t)) in 
f —)■ f -t- D{a{t), P{t)) with the prediction time dehned as 
cr(t) = t D[(j{t),P(t)). The key point of the design is the 


feasibility condition Tc dehned in ( |56l l, which ensures that the 
control action can reach the plant, namely, the delay rate is 
bounded by unity. The a priori satisfaction of ( [56| depends on 
the magnitude e and the angular frequency w of the periodic 
instability, and on the design parameters of the extruder. 

Theorem 3. For any setpoint x* £ (0, L) and for any choosen 
setpoint slope S{x*) £ K satisfying S{x*) > 5'min(a^*), where 
*S'min(a;*) is given by (113 and any initial condition xq £ [0, L) 
and, 

{Uo{e)\Uo{e)€[0,l), foralie&[-D{to,xo),0)}, (78) 

taking the control gains ai(x*) and ar(x*) as solutions of 
( |32p and (HI, respectively, the closed-loop system consisting 
of the plant (|77P-(|73P with state x{f), together with the 
control law ( |74P -( [77^ \22\ with actuator state U{t 9), 
9 £ [—D{t, x(t)),0), is GAS at x = x*, U = v{x*) if the 
parameters of the extruder model and the perturbation satisfy. 


or. 


(1 + 02 ^)' 


0 < 


< 


(l-e )2 


eu! 


< 


0102 


(1 + 02 ^) 2 ’ 


0 1 


(79) 


(80) 


0102 etu 40102 JO ^ 

(1 + 02L)2 ^ (l-e)2 < (1 + 02L)2’ 

( 81 ) 

where, 9i and 02 are defined in 0 and 0, respectively. 

Proof: The proof of Theorem is based on the Lyapunov- 
like condition ( [5^ that must be satished a priori to guarantee 
the GAS property for any given x* £ (0, L). In the following, 
we compute the function ( |49] l for the time- and state-dependent 
input delay model of the extruder (711-(73 1 , ( fTh) !, ( [TT] ) in order 
to establish that the feasibility region as it is dehned by ( |56l l is 
the entire physical domain, namely, x £ (0, L) and U £ [0,1). 
It holds that 

dD, . 0iea; sinfwf) , , 

= (82) 


Kt) 


dD 

dx 


{t,x)fit,x,U) =V{x,U), 


(83) 


where V{x,U) is given by 0. Since, 0 < e < 1, for all 
X £ (0,L) we get that minc(f) = 0i(l — e), and hence, 

dD . . euj(L — x) 

The gradient of ( [73| with respect to the input U satishes 

92 


SJuT{x,U) = - 


{l + 92x){l-Uf 


(85) 


It follows that ( [73| is a strictly decreasing function of U which 
belongs to \f),Vmax], for all x{t) £ [0,oo) and 

sup T{x,U) = . 

(1-f 02a:) 


( 86 ) 

































The delay rate ^ is uniformely bounded by unity, namely, 
the feasibility condition ([5^ is satisfied if and only if 


(87) 


< 1 . 


for all f > 0, X G (0, L) and U G [0,1). By (84i and (86i, 
it follows that (87i is satisfied for all x G (0, L) and U G 

[ 0 ; Xjfiax] 

1 ) if 

euj{L — x) 02X 


A(x) = 


(1 + 02 A) 

The derivative of the function 
A'(x) = 


< 1 . 


IS written as 
euj 


{1 + 92 x )^ 0i(l-e)2’ 


and 


(l-e)2 
sup A(x) = 


> 0102 , 


eujL 


a:G[0,L] 0l(l — 

Finally, the feasibility condition is satisfied if 

eui 


0102 < 


(l-e)2 L 


xi = (1 - e) 


01 


1 




and 


The relation 


A(xi) < 1. 
is equivalent to 
0102 


< 


euj 


< 0102 - 


relation 


ecu 


< 


(l_e)2 (1 + 02L)2' 


For satisfying ( |98| ) and ( |99l l we need to either impose 
or the condition ( |98l l with 02 < ^ which can be 
combined with (|9^ and (|94li in order to derive (|80|. ■ 


Remark 1. For given values for the parameters of the 


( 88 ) 


9i(l-e)2 (1 + 02X)’ 

is a strictly increasing function since then its maximum 
over the domain (0, L) satisfies 

02 ^ 


extruder, namely, 0i, 02, and L, the condition (79i is always 
satisfied if e or cu are sufficiently small. An increase of the 
magnitude of e causes a decrease of the allowed w and vice 
versa, as it is evident from (791- (HT). For given 0i, and L, the 
maximum bound of the perturbation parameters, namely, e or 
uj is expressed in ([79]l as 


sup 

e 2 eR 


010 - 


1172 


( 1 + 02^)2 


0 ^ 

4L’ 


02 = 


( 100 ) 


(89) 


(90) 


Larger variations of e and cu are possible, especially in the 
case in which 02 is small, as it is evident from ( |80l l and 
(811. However, for very large 02, one can conclude from 
( 7^ and iB that the allowable size of e and ui is restricted. 


and hence, since x G (0,L), (87i is guaranted if ( |79| ) 
holds. 

2) if ( [88l l is a decreasing function we deduce from ( |90l l that 
eui 


Moreover, from ([^-([IB one can conclude that the size of 
the allowable fluctuations of the transport speed in e and to is 
proportional to 0i and inversely proportional to the extruder 
length L. 


(91) 


(92) 


(93) 


In physical terms, conditions (79i- •HB are mainly a 
correlation between the pressure and the “rotation” flow. 


One should notice that ( |93| necessarily restricts 02 to 
satisfy 

02 < ^. (94) 

3) if ( [88] l is an increasing function of x on the interval 
[0,Xi] and a decreasing function on [xi,L] such that 
A(xi) < 1. The maximum value of A is attained at Xi 
satisfying 


(95) 


02 ew 02 

It becomes clear that ( [88] l admits a unique maximum 
satisfying the feasibility condition if 


namely, 02 defined in (151 and 0i defined in (14 1 , respectively. 
We recall the expression of the net flow rate defined in ([^ 
which is an increasing function of 02 as it is shown in 
Fig. 0 Therefore, changes in 02, by manipulating AT^, B, 
or po, the nozzle conductance, the screw resistance, and 
the melt density, respectively, affect the output flow rate 
Fd- For example, an increase in 02 by increasing the nozzle 
conductance Kd, leads to an increase in the outflow rate. Note 
that Kd, which defines the nozzle opening, is directly related 
to the printing resolution, namely, the accuracy of the printing 
process. A large nozzle opening leads to an extrusion of a 
filament with a large diameter and consequently deteriorates 
the printer precision. Moreover, from ( [80l l and ( [T5] l, it can 
be also seen that the “robustness” of the controller depends 
on the material thickness, namely, the mass density pg: a 
thicker material is less sensitive to large fluctuations of the 
transport speed under the predictor feedback control law. The 
parameter B in the expression of 02 in ( [T5] l is given by 

WH^ 

( 101 ) 


B = 


12 


(96) 

(97) 

(98) 


(1 + 02L)2 (l-e)2 

Using ( |88l l, the inequality ( |97] l leads to the following 

40102 

1 2 (99) 


where H is the approximate depth of screw channel from 
the screw thread root to the barrel internal surface and W is 
the width of screw channel. Consequently, changes in 02 due 
to the changes in B affect also the parameter 0i, since the 
screw pitch value ^ also depends directly on W. 

Relations (l^-lUB show that an increase in 0i, namely, 
an increase of the material convection speed, by enabling 
a large screw pitch ^ or a high screw speed iVg, improves 
the “robustness” of the controller in some way and allows 
for a system that supports broader changes of the convection 
velocity in frequency and amplitude. Note that a sharp 
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increase in the rotational screw speed Nq results in material 
overload and clogging problems and has a major effect on 
the residence time that is the critical time during which the 
material should be heated to have good properties before 
being evicted through the nozzle. Particulary, the extruded 
filament homogeineity is directly related to the residence time 
and to the process of solidification after layers deposition in 
3D printers. In addition, an increase in 9i in the screw speed 
Nq, increases the thermal energy in the extruder chamber due 
to the material shearing and decreases the viscosity of the 
melt. In that case, a rapid feeding of the extruder with granular 
material by applying a more agressive “Bang-Bang” control 
action absorbs the excess heat in the system. Maintaining a 
reasonnable temperature inside the barrel is essential because 
an excessive overheating of the system bums the polymer or 
produces poor extrusion. Generally, the conventional extrusion 
processes are equipped with a cooling system to compensate 
for the heat generated by the mechanical shearing effect that 
is proportional to the screw speed. 

In general, the nozzle and the screw designs are directly 
related to the predictor feedback control design and for 
achieving high performances for the closed-loop system the 
scale of the extruder should be neatly chosen. For instance, 
the agressiveness of the controller is influenced by the 
choice of 02 since the minimum value of the slope at the set 
point 5'inin(a;*) defined in •HI) depends on this parameter. 
Moreover, the entire process operates with an extruder head 
that moves very fast to print filament lines layer upon layer 
on a moving platform. A sufficiently light extruder head with 
small nozzle opening Kd and a small length L that operates 
at a sufficiently high screw speed Nq is needed to ensure a 
high rate of extrusion with a high precision. 


L = 0.2 m and the system is supposed to settle at Xq = 0.1m 
at the initial time. The value of the slope function is set to 
S{x*) = S'„iin(a^*) + 30. The simulations show the dynamics 
of the input filling ratio U{t), the interface position x{t), the 
predictor state P{t) and the delay function D{t,x{t)). Differ¬ 
ent cases including the open-loop dynamics, both uncompen¬ 
sated and compensated delay control laws are simulated for 
{e = 0.1, uj = 3.5 rad/s} and {e = 0.4, lu = 0.4 rad/s}. It is 
clear that, the uncompensated input leads to a limit cycle and 
the compensated closed-loop control allows faster convergence 
than the open loop control. Also, as it is shown in Figure]^ the 
feasibility condition is satisfied in both presented simulation 
results. 


S(i*) = 30, x’ = 0.16 , e = 0.1 ,u) = 3.5 



Fig. 5. Compensation of the time- and state-dependent input delay-(a). 



S(t*) = Smiu+ 30 , x* = 0.16 , e = 0.4 

w = 0.5 
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VII. Simulations 

A. Time- and state-dependent input delay compensation 

The setpoint is chosen as x* = 0.16 m which corresponds 
to a desired nozzle output flow rate as indicated by the 
equation (|^. The initial position of the moving interface 
is set to To = 0.1 m, the total length of the extruder is 




Fig. 6. Compensation of the time- and state-dependent input delay-(b). 


B. The state-dependent input delay compensator for the model 
with constant viscosity 


The following simulation results show the stabilization of 
the model described by (111 with the state-dependent input 
delay predictor feedback law | |20[ . Defining as 


x{t) 

Ds{x{t)) 


f {x,U {t - Ds{x{t))), 
L — x{f) 

0i ’ 


( 102 ) 

(103) 


where 


f{x{t),U{t)) 


e^T{x{t),u{t)), 


(104) 
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9i is the nominal transport velocity of the material defined 


in (14i, and the function T(x{t)^U{t)) is given by ( [7^ . The 
predictor feedback controller is written as 


/(P,(p),C/(m)) 


(105) 
dp. (106) 


P,{t) = x{t) + [ T _ p (p / \ rrf \\ 

where for all t — Ds{x{t)) < fJ- <t 
QD 

Fs{Ps{^i),U{^l)) = ^(P,(p))/(P,(p),[/(p)). (107) 

By specializing Theorem to the case e = 0 it can be 
shown that the predictor feedback law (105 i- (107 1 renders 


S{x-] = 30 

x' = 0.16 



- - Open-loop 
■■•Uncompensated 



—- 

.Predklor 


<r(0)y^[P(0)) ... 



2 4 

6 8 10 12 1 

4 16 

18 20 


More precisely, both the predictor state, P{t) and the delay 
function, D{x{t)), are estimates of the actual prediction state, 
namely, P{t), and delay function, namely, D{t,x{t)), that 
are described by an d ([T6 ] i, CD’ respectively. For im¬ 
plementing the controller ( |108| ), ( |109|l, t he “actual” feasibility 
condition, defined in (|79]l, ( |80| ), and (|81[l have to hold, in order 
to guarantee that the controller actually “kicks in”. In addition, 
we assume that the following condition, which guarantees that 
the denominator in ( |109| l remains always positive (and hence, 
the controller remains bounded) is satisfied 


e2Pio) 


u{e) 


system ( |102| l-( [T04l l GAS (in the physical domain) at any given 
equilibrium x*. 


(l+ 02 P( 0 ))(l-C/( 0 )) 
0 e [f — D{x{t)), t]. 


(1-C7(0)) 


< 1 , 


(113) 


Note that with strictly positive physical parameters B, pQ, and 
Kd, and for ?7 G [0, Vmax], relation (113 i is satisfied whenever 


P(0) € [0, oo), for all 9 € [t — D{x(t)), t]. 

The simulation results in Fig. illustrate that the state- 
dependent input delay compensator can handle small time- 
varying uncertainties on the vector field ( | 112 [ ) and the delay 
function ( |lll| i, as described by ( 1 ^-(|^ and ([T§-([T7|- 
respectively. An increase in the necessary control effort to 
drive the system to the setpoint is also denoted and the rate 
of convergence decreases compared to the time- and state- 
dependent predictor feedback (|7^-(|76ll shown in Fig. 


S{x’) = 30 , x’ = 0.16 , e = 0.4 ,L3 = 0.5 


Fig. 7. Compensation of the state-dependent input delay. 


C. Control with a state-dependent input delay compensator 

We deal with the case in which the time variations of 
the transport speed are unknown and consider the closed- 
loop system consisting of the plant ( [T3] l with an actual delay 
D{t, x{t)), given in (16 1 , (ID, together with a state- dependent 
input delay predictor feedback defined by 


Uit) = V (P{t)) 
P{t) = x(t) + 


f{Pis),Uis)) 

t-b(x(t)) l-P(P(s),P(s)) 


(108) 
ds, (109) 


where. 


dD 


P(P(s),P(s)) = _(P(s))/(P(s),P(s)). 

With an estimated delay function defined as 

L — x{t) 


D{x{t)) = 

and the nominal vector field, 

f{x{t),U{t))=9, 

-f 


^No ’ 
92x{t) 


( 110 ) 


( 111 ) 


{l + 92xit))il-U{t)) 
Uit) 



18 20 


Fig. 8. Robustness of the state-dependent input delay compensator. 


5(x-) = 



1-Uit) 


U(t) G [0 1) (112) Feasibility condition: the delay rate is bounded by unity. 
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VIII. Conclusions 

This paper is devoted to the stabilization of a screw- 
extrusion process. A coupled PDE-ODE model is used to 
derive a state-dependent input delay system describing the 
melt convection in the extruder chamber for an isothermal 
case. The extension of the aforementioned model to a non- 
isothermal case is proposed introducing a periodic time- 
dependent function in the state-dependent input delay function. 
Next, we design a predictor feedback controller to compensate 
the state- and time-dependent input delay and establish the 
GAS of any setpoint with respect to the physical domain 
under physical and design restrictions. The delay compensator 
is constructed with a nominal Bang-Bang-like controller that 
ensures the GES of the delay-free plant. 

It is clear that the delay function model for the non- 
isothermal extrusion process should depend on the rheological 
properties of the extruded polymer. In general, a fairly 
accurate estimation of the material friction and the viscosity 
behavior is extremely hard to achieve in such processes 
due to the change in material composition, and the strong 
interaction between the heat and mass transfer phenomena. 
An interesting future work would be to consider an unknown 
time-dependent perturbation acting on the polymer convection 
speed. As it is shown in the simulation results, the state- 
dependent compensator might be able to handle time-varying 
perturbations acting on the vector field and the delay function. 
The proof of the robustness properties of the state-dependent 
input delay predictor is an interesting perspective. 

Physical definition of the parameters 


system 

i{t) = f + W{t- D{t,x{t)))) (A-2) 

Wit) = 0. (A-3) 

Proof: The proof of Lemma [T] is based on a di¬ 
rect verification considering P {t — Dit,x{t))) = x{t) and 


a {t — Dit, x{t))) = f in the original system (431. 


Lemma 2. (Inverse Backstepping Transformation) The in¬ 
verse of the infinite dimensional backstepping transormation 


[A-11 is defined for all t — 0(1, x(t)) < 9 <t by 
UiO) = Wi9) + ^idi9),Uie)), 


(A-4) 


with 


n(0) = 


/((T(s),n(s),K(tT(s),n(s)) -f lE(s)) 


/0(t) 

+ x(t) 

F(a(s),n(s),VE(s)) 
do 


l-F(a(s),n(s),lE(s)) 


do 


^ Ik (c^(s),n(s)) 


X / (o-(s),n(s), /c((t(s), n(s)) 
1 


r(9) — t-i-f 
J d) 


l-F(a(s),n(s),tV(s)) 

(j){t) = t — D(t, x(t)). 


lV(s)) 

ds, 


ds 

(A-5) 

(A-6) 


(A-7) 

(A-8) 

(A-9) 


L = 200.10-3 

m 

Extruder length 

li 

O 

rpm 

Screw speed 

B = 9.3450 10-9 

4 

m 

Geometric parameter ^ 

Fd = — 

Kg s-i 

Net forward mass flow rate 

Kd = 2.45 10-3 

m? 

Nozzle conductance 

^ = 10 10-3 

m 

Screw Pitch 

Seff = - 

mf 

Effective area 

^eff = ^Seff 

rnf 

Effective volume 

v = — 

Pa s-^ 

Melt viscosity 

Po = 1240 

Kg m-3 

PEA Melt density 

e = - 

— 

Amplitude of the perturbation 

UJ = - 

rad.s~^ 

Erequency 


Proof: Direct verification considering that P(6) = 11(0) 
and d(9) — a(9) for all t — D(t,x{t)) < 9 < t. We refer 
to P(9) as the plant-predictor system and 11(0) as the target- 
predictor system, respectively. However, they play different 
roles because they are driven by different inputs (U versus 
W). 

Lemma 3. (Stability of the Target System) For any positive 
constant g, there exist a class /Coo function p* and a class ICC 


the feasibility condition the following holds: 

m < 

m = 


/3(p*(-(/o)),/-/o), t>to 
|a;(f)|-I- sup |VE(0)| 


where 


p*(s) = 




(A-10) 

(A-11) 


(A-12) 


APPENDIX 


We recall Lemmas 1-8 from p()[ which are applied to the 
nonlinear time- and state-dependent input delay system ( |43| ) 
as Lemmas [gig for the proof of Theorem 

Lemma 1. (Backstepping Transformation of the Actuator 
State) The infinite dimensional backstepping transform of the 
actuator state given by 


1 — c 

Proof: Based on the input-to-state stability of i; = 
/ (t, X, K(t, x) -f ijj) with respect to w, namely, Assumptionj^, 
there exist a smooth function S(t, x{t)) : [to, oo) x K" —)■ IR_|_ 
and class /Coo functions ai, a 2 , as, a 4 such that 

ai(|a;(/)|) < S{t,x(t)) < a 2 i\x(t)\), (A-13) 

S(t,x(t)) < -a 3 (\xit)\) aiilWit - D{t,x(t)))\), 

(A-14) 


with 


1V(0) = C/(0)-«((t(0),P(0)), 


(A-1) Sit,xit)) = 


for all t — D(t, x(t)) < 0 < allows to transform the system 
with the controller into the following target 


dS{t,x(t)) dS(t,x(t)) 
dt dx 

X / (t, x(t), K(t, x{t)) -CWit- D(t, x{t)))). 

(A-15) 
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Let us define the Lyapunov function for the target system 0 
and ( |A-31 l as 


where 

m = 


pL{t) / \ 

V(t) = S(t, x(t)) + k / - dr, (A-16) 

Jo ^ 


sup 

t — D{t,x{t))<d<t 

= lim ( / e 

n->oo V Jt-D{t,x{t)) 


g{l+a{e)- 




2ng{l-\-(T{9) — 




l/2n 


(A-17) 


with g > Q. Let us upperbound and lowerbound ( |A-17| l in 
terms of: 


sup 

t-D{t,x{t))<e<t 


Wid) 


(A-18) 


Upperbound of L(t): Using the feasibility condition (56 1 
and (|55|), we deduce 


m < ^ 

i — I 


(A-19) 


By integration of (A-191 on [t — D(t,x{t)) 9] with 
a {t — D{t,x{t))) = t, we derive the inequality 

l + cr(6 »)-f < j^{l +D{t,x{t))), 

y t-D{t,x(t)) <9 <t. (A-20) 

From Assumption and (j 6 ^, the following inequality 
holds: 

L{t) < sup \W{9)\. 

t — D(t,x{t))<9<t 

(A-21) 

Lowerbound of L{t): Similarly, using the fact that 
a {t — D{t,x{t))) = t , with 9{t) being an increasing 
function, we obtain 


l + a{0)-t > 1 , 

and hence 

L{t) > e® 

V t —D{t,x{t)) <9 <t. 

The time derivative of (|A-17[l is 


sup \W{9)\ 

t—D{t^x{t))<9<t 


(A-22) 


(A-23) 


L{t) = lim — 
n—>oo 2n 


t l/2n—l 

2r^g(l+^T(fl)-t)TJA/'a^2n,- - ' 




>W{9Y^d9 


I - (1 - - D{t,x{t))f^ 

I at 

2ng [ (A-24) 


t - 

< 1 and hence 


where = t — D{t,x{t)). By (56 1 , it is clear that 

dD{t^x{t)) 
dt 


L[t) < -gL{t). 


(A-25) 


Computing the derivative of the Lyapunov function ( A-16 1 as 


V{t)= S{t,x(t)) + kL{t)- 


m ’ 


we deduce 

V(t) < S{t, x(t)) — kga 4 {L(t)). (A-27) 

Using the boundness of S{x{t)), ( |A-14 i, the following in¬ 
equality holds: 

V{t) < -a 3 {\x{t)\) + a 4 {\W{t - D{t,x{t)))\) - kga 4 {L{t)). 

(A-28) 


Imposing k = g by (A-23 i, we derive the inequality 

V{t) < -a3{\x{t)\) - a^mt)), (A-29) 


and with ( |A-13| l, ( |A-16| l and ( |A-17| l, we conclude that there 
exists a 1C function 71 such that 


m < 


(A-30) 


By the comparison principle, there exists a class JCL function 
/3 such that 

V{t)<P(Vito),t-to). (A-31) 

From ( A-13[ ) and ( |A-16| l and the properties of class JCL 
functions, we finally get 

l^(^)l + ^ IJi{\x{to)\ + — to), (A-32) 

/3i being a class ICC function. Considering the boundness of 


L{t) defined in (A-20i and (A-23i, the prove is achieved. 


Lemma 4. (Bound of the Predictor in Terms of Actuator 
State) There exists a class ICao function p such that for all 
the solutions of the system satisfying the feasibility condition 


(561, the following holds for all t — D{t,x{t)) < 9 <t 


|f’( 6 ')|<P |a:(f)|+ sup |C/(s)| • (A-33) 

Y t—D{t,x{t))<s<t J 

Proof: Differentiating ( |47] ), we deduce the following rela¬ 
tion for all t — D(t, x{t)) < 9 <t 

dP{9) ^ f{a{9),P{9),U{9)) 

d9 I- F(a{9),P(9),U{9)y 

F(a{9),P(9),U{9)) = {a{9), P{9)) f{a{9),P(9),U{9)) 


dD 

-^{<J{9),P(9)) 


(A-34) 


and with the change of variable y = a(9), (A-34i may be 
rewritten as: 

= f (y, P( 0 (y)), U{y - D(<P{y)))), (A-35) 

ay 

t<y< a(t). 

From Assumption [T] we get that 

dR{y,P{(l){y))) 


d9 


<&{9) {R{y,P(^(y))) 

+ P 3 (\U{y-D{ck{y)))\)), (A-36) 


for all t < j/ < a{t) and using the feasibility condition (56i, 
we deduce, for all t — D{t, x{t)) < 9 < t 


d9 


1 — c 


(A-37) 


(A-26) 
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By Assumption and the comparison principle, we obtain 

R{a{e),P{e)) < 

(R{t,x{t))+ sup ^ 3 (|C/(s)|) ] , 

\ t — D{t,x{t))<s<t j 

t-D{t,x{t)) <e <t. (A-38) 


With the standard properties of class /Coo functions the Lemma 
0 is deduced and the class ACoo function p is written as; 

p{s) = /rr' ((M 2 (s) + ri3(s)) . (A-39) 


Lemma 5. (Bound of the Predictor in Terms of Transformed 
Actuator State) There exists a class K, function ^ such that 
for all the solutions of the system satisfying the feasibility 
condition the following holds: 

|n(6»)| < i/i ( |a;(f)| + sup |VF(s)|], (A-40) 

Y t—D{t^x(t))<s<t J 

for all t — D{t, x(t)) < 9 < t 

Proof: The plant X{t) = / (t, x{t), K(t, x{t)) + Lu(t)) satis¬ 
fying the uniform input-to-state stability property with respect 
to w, and the function k being locally Lipschitz in both 
arguments and uniformly bounded with respect to its hrst 
argument, there exist a class ICC function /32 and a class AC 
function -fi such that for all t > to 

Y{t)< ^2(|T"(fo)|,'r-fo) + '0i fsup |u;(s)|V (A-41) 

\s>to / 

with 


Y{t)= f (Y{t),k{t,Y{t))+uj(t)) . (A-42) 

Now, we consider the change of variable y = a{9) and write 
the predictor of the target system (A-5 i as 

du(m) 


dy 

t<y< o(t). 


= f iy, n((?!)(2/)), k(?/, n((;i(i/))) + w{(/)(2/))), 

(A-43) 


where Cl and 2 are defined in ( |64| l and ( | A-11 1 ), respectively. 

Proof: Using the inverse transformation ( A-4| l and the 
bound ( |A-40| l, we derive ( |A-45| l with 


di (s) = s + /5(V'(s)), 


(A-47) 


and from the direct transformation ( |A-l| i together with the 
bound (|A-33|l, we deduce (|A-46[), where pi is dehne as 


px{.s) = s + p[p{s)) (A-48) 


Lemma 7. (Ball Around the Origin Within the Feasibility 
Region) There exists a positive constant 7 such that for all 
the solutions of the system that satisfy 

\x(t)\+ sup |17(6I)|<7, (A-49) 

t — D{t,x{t))<9<t 

the feasibility condition is satisfied. 

Proof: From ( |45] l we derive the following inequality 

\f {t,x(^t),U(t - D{t,x(t))) I < 

a ^|x(f)| + supi_£,(t_^(i)) 

<s<i |C/(s)|) (A-50) 

Recalling the relations and ( [62] i of Assumption we 
deduce that for all 9 G [t — D(t,x{t)),t] and c €]0,1[ , if a 
solution satishes 

c, + pfy\P(9)\) + {c 2 +pe{\m\)) 

Q! ( |P( 6 >)|-f sup |C^(s)|) < C, (A-51) 
Y t — D{t,x{t))<s<t J 

then it also satishes 

Using Lemma we conclude that ( |A-51| l is satished if the 
following holds 

(c 2 + peipi^it)))) a. (p(U(f)) -f Cl(t)) 

+ p 5 {p(Cl(t))) < C-C 3 . (A-52) 

Let us dehne a class ACoo function pc as 

Pc{s) = P5{p{s)) + (C 3 + P6(p{s))) a {p(s) + s) . (A-53) 


Using ( A-42| l, we derive the following relation 

|n( 6 »)| < fi2{\x{t)\) + 


sup |VF(s)| , 

\t—D{t,x{t))<s<t J 

(A-44) 

for all t — D{t,x(t)) < 9 <t with a class AC function ^’ 2 ( 5 ) = 
fi 2 {s, 0) . Using the properties of class AC functions, (A-40l is 
deduced with tpis) = tpiis) + ' 02 (s)- ■ 


Lemma 6. (Equivalence of the Norms of the Original and 
the target system) 

There exist class ACoo functions pi, pf such that for all the 
solutions of the system satisfying the feasibility condition 
and for all t > to, the following hold: 

n(t) < pj\E(t)), (A-45) 

2(t) < pi{n{t)), (A-46) 


It follows that 

7 =/O0^(c-C 3 ). (A-54) 


Lemma 8. (Estimate of the Region of Attraction) There exists 
a class 1C function '0roA such that for all initial conditions of 
the closed-loop system that satisfy relation (63 i, the solutions 
of the system satisfy (A-49 \ for c g]0, 1[, and hence, satisfy 

Proof: Using Lemma and (|A-10[), the following holds: 


Cl(t)<pT (/3(p*(pi(U(fo))),A-Ao)), 


(A-55) 


where U is dehned in ( |64| ). Introducing the class ACoo function 
pg{s) = P 7 ^ (/3(s,0)), we derive the inequality 


^(A) < P9 {p*{pi{^{to)))) ■ 


(A-56) 
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Hence, for all initial conditions that satisfy the bound ( |5^ 
with any class 1C choice 

ipRoAic-ca) < Pi'^{p7HtJ'9^iPc\c- C 3 )))), (A-57) 

the solutions satisfy ( |A-49| l. Moreover, for all of those ini¬ 
tial conditions, the solutions verify ( [64| ), for all 9 > to — 
D{to,x{to)). m 
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